Full Counting Statistics for Orbital-Degenerate Impurity Anderson Model 
with Hund's Rule Exchange Coupling 
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We study nonequilibrium current fluctuations through a quantum dot, which includes a ferro- 
magnetic Hund's rule coupling J, in the low-energy Fermi liquid regime using the renormalized 
perturbation theory. The resulting cumulant for the current distribution in the particle-hole sym- 
metric case shows that spin-triplet and spin-singlet pairs of quasiparticles are formed in the current 
due to the Hund's rule coupling and these pairs enhance the current fluctuations. In the fully 
screened higher-spin Kondo limit, the Fano factor takes a value Ft, = (9M + 6)/(5M + 4) deter- 
mined by the orbital degeneracy M. We also investigate the crossover between the small and large 
J limits in the two-orbital case M = 2, using the numerical renormalization group approach. 

PACS numbers: 71.10.Ay, 71.27.+a, 72.15.Qm 
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The experimental study of electron transport in meso- 
scopic devices, subject to applied bias voltages, has pro- 
vided a new way of investigating interelectron many-body 
effects pHil . For example, the theoretical prediction that 
the Kondo correlation on the nonequilibrium currents 
through quantum dots leads to an enhancement of the 
shot noise in these systems 0-0] has stimulated subse- 
quent experiments [8l4l0l|. 

A deeper understanding of nonequilibrium transport in 
these systems can be obtained from the calculation of the 
current distribution function, the higher order cumulants 
beyond averaged current, and the noise power. These 
can be calculated from the cumulant generating func- 
tion (CGF) for nonequilibrium transport. This quantity, 
however, is difficult to calculate when many-body effects 
have to be taken into account. A recent development 
has been the calculation of the current probability dis- 
tribution for quantum dot in the Kondo regime 11-13} 
using the general formulation of the full counting statis- 
tics (FCS) [3, [HI ■ It was shown that quasiparticle sin- 
glet pairs carry a charge 2e in the backscattering current 
due to the Kondo correlation and that these enhance the 
nonequilibrium current fluctuations. 

In this Letter, we investigate the FCS for an orbital- 
degenerate impurity Anderson model as a prototype 
model to examine the effects of both the Hund's rule 
and Kondo correlations and, hence, to deduce the time- 
averaged current and shot noise. These have been experi- 
mentally investigated in vertical dots, carbon nanotubes, 
and double dots 0, [1(3, H3] • Using the renormalized 
perturbation theory (RPT) [lij], we calculate the zero- 
temperature CGF up to third order in applied bias volt- 
age for arbitrary strength of the interactions at the dot 
site in the low-energy Fermi-liquid regime. We show that 
the Hund's rule correlation gives rise to spin-triplet pairs 
and spin-singlet pairs of quasiparticles carrying charge 
2e in the nonequilibrium current, which characterize the 



shot noise and higher order cumulants at low energies. 

Model. — Let us consider transport through a corre- 
lated dot described by an orbital-degenerate impurity 
Anderson model H A = rl + H T + Hj with 

%0 = £ ka c kamcr C kama- + ^ e d n dm<j > (1) 

karaa ma 

+ H.c), (2) 

kama 
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Here, d ma annihilates an electron in the dot level e d 
with orbital m — 1,2,- • • ,M and spin a, and c kam<J 
annihilates a conduction electron with moment k, or- 
bital m, and spin a in lead a — L , R. Interactions U, 
W(0 < W < U), and J(< 0) are the intra- and in- 
terorbital Coulomb repulsion, and Hund's rule coupling, 
respectively. n dmiy is the number of electron in the dot 
level with orbital m and spin a, and S dm is the spin 
operator of electrons in orbital m of the dot. The in- 
trinsic level width of the dot levels owing to tunnel cou- 
pling v a is given by T = ^ Q i^p c v 2 a with the density of 
state of the conduction electrons p c . The two-orbital case 
(AI = 2) has been experimentally realized in such sys- 
tems as vertical dots, carbon-nanotube dots, and four- 



terminal double-dots [l6, 17 1. For simplicity, the sym- 
metric lead-dot coupling v L — v R and the particle-hole 
symmetry e d = —U/2 — (M — 1)W are assumed. The 
chemical potentials of the leads P L / R — ±V/2, satisfying 
fi L — [i R — V(> 0), are measured relative to the Fermi 
level defined at zero voltage V = 0. We take units with 
h = fee = e = 1 throughout this Letter. 

Full counting statistics. — The probability distribution 
P(q) of the transferred charge q across the dot during a 
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time interval T can provide the current correlation func- 
tions to all orders. In order to treat them systemati- 
cally, we calculate the CGF In* (A) = ln^, e lXq P{q) 
in the Keldysh formulation 0: In* (A) = hi (T C S&) 
where S£ = T c exp{-i j c dt [H^(t) + Kj(t)] } is the 
time evolution operator for an extended Hamiltonian 
= Hq + Hx + Hi, C is the Keldysh contour along 
[t : -T/2 -> +T /2 -> -T /2], T c is the contour-ordering 
operator, and A is the counting field. Here, is given 
by 



w* = E 

kma 



e lX(t)/2 d] na c kLma + iVLc fciw + H.c(k) 



with the contour-dependent counting-field defined by 
A(i) = = ±A for the forward and backward paths 
labeled by "— " and "+" respectively. 

To calculate the CGF, we use Komnik and Gogolin's 
procedure [2(| outlined below. First, a more general func- 
tion %(A_,A + ) is introduced. It is basically given by 
In (TcSq) but A T is formally treated as an independent 
variable assigned for each contour. For the long time 
limit T — > oo where the switching effect is negligible, 
the general CGF is proportional to T. Then, the deriva- 
tive of the CGF with respect to A_ is given in terms of 
Green's functions as 

<lnx(A_,A + ) 



krna 

with A = A_ — A + . .g°« + M = i27T(5(w - E ka )f a {u) 
and = -i27r5(cj - e fca )[l - are the lesser 

and greater parts of the Green's function for electrons in 
lead a, respectively, with the Fermi distribution function 
f a (w) = [e^ u ~i J -o,)/ T + For the long time limit T -> 
oo, the dot Green's function is defined as G\ vv (uS) — 

-^^-Oe-C-O^C^jiftOlA^v " and 
v 1 are the labels for the two Keldysh contours. Here, 
we use a notation (A(t)) x jA = (T c S&A(t)) /x(A_, A + ) 

which represents an expectation for Hamiltonian T~L\- 
Once lnx(A_,A + ) is computed, the statistics are recov- 
ered from lnx(A) = lnx(A, —A). 

Renormalized perturbation theory. — To calculate the 
Keldysh Green's function with finite A at low ener- 
gies, we use the RPT outlined below [H SH, [H| . The 
basic parameters that specify the impurity Anderson 
model W A are e d , T, U, W, and J. Correspondingly, 
the low-energy excitations can be characterized by the 
renormalized parameters for the quasiparticles: dot-level 
e d = z[e d + TT d (0)], level width f = zT, where Y> r d {u) 
is the self-energy of the retarded Green's function for 



[1 — dYi r d (uj) /du}\ u= o] 1 is the wave function renormaliza- 
tion factor. Here, T corresponds to the characteristic en- 
ergy scale, namely, the Kondo temperature Tk = 7rF/4. 
The quasiparticle interactions U, W and J are defined 
in terms of the local full four-vertex for the scattering 
of the electrons, T^IZV.^ ^H^ , ujy, uj^, U4) taken at zero 
frequency u>i = (see Ref. l2lh. In this Letter, we choose 
the renormalized parameters defined at the equilibrium 
ground state. Then, the replacement of the bare param- 
eters with the renormalized ones gives the leading terms 
of an effective Hamiltonian corresponding to the low- 
energy fixed point of the numerical renormalization group 
(NRG). The renormalized parameters can be related to 
the enhancement factors of spin, orbital and charge sus- 



ceptibilities; zx s = 1 + p d (0) U - (M - 1) J 



2W-U 



,ZXc = l ~Pd(°) 



of 



U + 2(M- 1)W 
state Prf(w) 



1 + ^(0) 

with the renormalized density 
(T/tt)/[(uj — <f d ) 2 + T 2 ], and also to the Friedel sum rule 

7r(n dmcr ) = cot -1 (e d /Tj. These relations enable one to 
deduce the value of renormalized parameters in some spe- 
cial limits. For instance, the renormalized dot level in 
the particle-hole symmetric case is situated at the Fermi 
energy 1 d = 0. We can also evaluate the renormalized 
parameters for dots with two orbitals (M = 2) with the 
NRG approach [p. 

The nonequilibrium perturbation theory in powers of 
U, W and J can be reorganized as an expansion with 
respect to the renormalized interactions U , W, and J by 
taking the free quasiparticle Green's function of the form 



H = 
(«) = 



'V, 



iTe fx ' 2 f L + iTf h 



9 d 



*+-(«) = 



(oj) = -oj + e d -iT + iT J2i ft 



'V, 



(0) 



with D 



( w -? rf ) 2 + r 2 + r 2 ( e tf/a_i) (i-f R )f L 



the zero-order propagator [22j. With this approach, the 
exact form of the Green's function can be calculated at 
low energies up to terms of order u 2 , V 2 , and T 2 . 

We calculate the Green's function G d (uj) under a finite 
counting field through the Dyson equation 



(7) 



The renormalized self-energy at T — up to ui 2 , uV , and 
V 2 in the particle-hole symmetric case can be calculated 
in the second order perturbation in the three renormal- 
ized interactions as 



the dot state G r d {uS) 



+ W - TT d ( U )]-\ z 



sr 



A(u,V) B(u,V) 
-B*(-u,-V) A(u,V) 



(8) 
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with 



A{uj,V) 
B(u,V) 



a(ui,3V/2) + 3a(w,F/2) 

+ 3a(w,-V/2) + a(w,-3V/2)], (9) 

e~ a ' 2 b (w, 3V/2) + 3& (w, V/2) 

+3e iS / 2 6 (w, -V/2) + e iS 6 (u, -3V/2)1(J0) 



S 2 



+ 2(M-1) (w 2 +3j 2 /4) 



(11) 



Here, a(w, x) = — j(w — x) 2 sgn(— uj + x), b(tu, x) = (cj — 
:r) 2 0(-u; + a:), and it ee U/(nf),w = W/{ttT) and j ee 
J/(tiT). For A = 0, it corresponds to an extension of Ref. 
to the orbital-degenerate Anderson model. 

Results and discussion. — Integrating Eq. ([3]) with re- 
spect to A_, the CGF is calculated at T — up to V 3 , 
as lnx(A) = F Q + F l + F 2 with 



2MT 



2tt 

Tx = TP b i 



V/2 

-V/2 
-iX 



dujln [l + T{u) (e iA - l)] , (12) 



1) , Ti = TP b2 (■ 



-i2A 



1) -(13) 



Here, J- Q describes free-quasiparticle tunneling via the 
renormalized level with transmission probability T(uS) — 
r 2 /(ui 2 + T 2 ), which carries single charge e. F 1 and T 2 
represent the^scattering process due to the residual inter- 
actions U, W, and J. From the coefficient of the count- 
ing field in F 1 and F 2 , it is clear that Pb\ = j^l^- 

and Pb2 = are the probabilities of the single- 

and paired-quasiparticle backscattering processes carry- 
ing charge e and 2e, respectively. Then, the nth-order 
cumulant C n = (— i) n ^pr hax(A) is readily derived as 

C n = T [Um + (-l)"(Ao + Pbi + 2 n P b2 )} , (14) 

where I u = 2MV/(2ir) is the linear-response current and 
P&O = represents the probability of the single- 

quasiparticle backscattering processes due to the renor- 
malized level. Especially, the factor 2™ characterizing 
cumulant of interacting electrons (fT4l indicates the ex- 
istence of spin-singlet, orbital-singlet, and spin-triplet 
pairs of quasiparticles in the nonequilibrium current. 
Equation (|14[) generalizes the previous results of the cu- 
mulant for the SU(2M) case where U = W and J = 
11- Til]. In the present result, the Hund's coupling enters 
through the P b \ and P b 2 as well as the Kondo energy scale 
F. Therefore, it gives rise to the spin-triplet pairs in the 
current, and also varies the couplings of the spin sector u 
and the orbital sector w, which causes the crossover seen 
in the current fluctuation as discussed below. 

In the series of cumulant C n , the first two coefficients 
C Y and C 2 give time- averaged current I = C 1 /T and shot 
noise S = 2C 2 /T, respectively. The higher order ones, 
C n for n > 3, are determined by these two and are auto- 
matically generated by Eq. (|14l) . From these coefficients, 



the Fano factor of the shot noise for the backscattering 
current ee I u — I, can be deduced in the form 



F s 



1 + 91 
1 + 51' 



(15) 



This is an important result of this Letter [23J . The Fano 
factor Ff, determines the effective charges of the backscat- 
tering current [4]. Our results, Eq. (fT5|) with (fTTj). show 
that F b also represents the degree of the residual inter- 
action, as the Wilson ratio R — z\ s does in the SU(2) 
case. We examine the value of the Fano factor in some 
special limits in the following. 

First, in the weak coupling limit it, w, \j\ <C 1, the 
renormalized parameters take a form (u, w, j) — > (u, w, j) 
and z -> 1 - (3 - \tt 2 ) [u 2 + 2(M - 1) (w 2 + 3j 2 /4)] 
with scaled interactions u ee U/(ttT),w = W/(ttT), and 
j ee J/(nT). Substituting these into Eq. (fl5)) . the re- 
sult for the second-order perturbation in the bare inter- 
actions is produced. In particular, the Fano factor for 
noninteracting limit u = w = j = naturally becomes 
unity which represents single charge transport given by 
Eq. ([121). 

In the large Hund's coupling limit \j\ ^> 1, the dot 
state is restricted to the highest spin state, for which we 
obtain the Kondo coupling H K from the Anderson model 
T-La in the form 



H 



K 



„t 



s d 



(16) 



kk' 



with J K = 4(vl + v 2 R ) [U - (M - 1) J] , the Pauli ma- 
trix a aa >, c kma = Y. a c k am J^ and the S = M/2 spin 
operator for the dot-site Sd- In this limit, the high spin 
S = M/2 (> 1) state is fully screened, and the Wilson 
ratio becomes R = 2{M + 2)/3 [2J]. The charge and or- 
bital fluctuations are suppressed Xc ~+ and Xorb — ^ 
in the limit j — > — oo. Therefore, the renormalized pa- 
rameters take the value (u,w,j) — > (1,0,-2/3). Then, 
through Eqs. (fTTj) and (fl"5j) . we find the Fano factor in 
this limit, 



9M + 6 _ 9S + 3 
5M + 4 ~ 5S + 2 



(17) 



We note that to observe the Fano factor due to the 
higher-spin Kondo effect in real experiments, Hund's cou- 
pling J is not required to be very large as long as it is 
larger than the Kondo temperature J > Tk- The ex- 
plicit values of the Fano factor in this limit for several 
M are shown in Table HI For comparison, the Fano fac- 
tor in the SU(2M) Kondo limit u = w — > oo,j = 0, 
F b -> (M + 4)/(M + 2) @,H3, is also shown in Table ffl 
Increasing orbital degeneracy M in the SU(2M) Kondo 
limit, the system approaches the mean-field limit Fb — > 1 
13]. In contrast, for large the renormalized interac- 
tions u, w, and j converge to values independent of the 
orbital degeneracy M. As M increases, the number of 



4 



TABLE I. The Fano factor F b = C 2 /(TI b ) in (a) the large 
Hund's coupling j — > — oo or the S = M/2 Kondo limit, and 
(b) the SU(2M) Kondo limit u = w -¥ oo, j = 0. 
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— > oo 


(a) S = M/2 Kondo 


12/7 


33/19 


7/4 


-> 9/5 


(b) SU(2M) Kondo 5/3 


3/2 


7/5 


4/3 


-> 1 




FIG. 1. (Color online) (a) The Fano factor for several choices 
of Coulomb repulsion u = w as a function of scaled Hund's 
coupling j. (b) The renormalized parameters for u = w = 3 
and (c) u = w = 0, as a function of j. 



coupled orbitals with the Hund's coupling J increases. 
It results in an enhancement of the spin fluctuations, as 
seen in the expression of z\ s or the interaction factor 
of the self-energy given in Eq. (fTTj). Therefore, we con- 
clude that the Fano factor increases with M, through the 
enhancement of the quasiparticle-pair scatterings due to 
the residual Hund's coupling j. 

Next, J-dependent crossover of the Fano factor for the 
two-orbital case (M = 2) is investigated, evaluating the 
renormalized parameters with the NRG approach. The 
Fano factor for u — w and u — 3 > w as a function of 
Hund's coupling j are shown in Fig. QJa) and Fig. EJa), 
respectively. The crossover is observed in the Fano factor 
around the coupling corresponding to the Kondo temper- 
ature |j| ~ tjc = r/(4r). At large \j\, the Fano factor 
for any u > w converges to a value 12/7 given by Eq. 
(fT"7|) . as discussed above. However, at small \j\, the val- 
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FIG. 2. (Color online) (a) The Fano factor for several 
Coulomb interaction w(< u = 3) as a function of scaled 
Hund's coupling j. (b) The renormalized parameters for 
(u, w) = (3, 0) as a function of j. 
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FIG. 3. (Color online) (a)The Fano factor and (b) the renor- 
malized level width for no Hund's coupling j — and several 
intraorbital Coulomb repulsions u, as a function of interor- 
bital Coulomb repulsions w(< u). 



ues of the Fano factor depend on Coulomb repulsions u 
and w. The Fano factor and the renormalized level- width 
(or the Kondo temperature) at j — for several intraor- 
bital Coulomb interaction u as a function of interorbital 
Coulomb interaction w are shown in Fig. [3] (a) and (b), 
respectively. With an increase of interorbital Coulomb 
repulsion w, the Fano factor varies from the SU(2) value 
to the SU(4) value around (u — w) ~ tx- 

The crossover seen in the Fano factor is a consequence 
of the behavior of the renormalized parameters, which 
enter through Eq. (IT5|) . Thus, the feature discussed above 
is also seen in the renormalized parameters, which for 
(u,w) = (3,3), (0,0), and (3,0) are shown in Figs. [T](b) 
and (c), and Fig. [2] (b), respectively, as a function of j. 
Particularly, in Fig. []] (b), we see a general trend that 
the spin coupling of the Kondo correlation u increases, 
whereas the orbital coupling w decreases, with increase 
of the ferromagnetic coupling — j. These renormalized 
parameters converge to the value ({t, w, j) — (1, 0, —2/3) 
in the large Hund's coupling limit. Even for (u, w) = 
(0, 0), the Hund's coupling enhances the spin coupling u 
as shown in Fig. QJc). 

Summary. — We have studied the role of the ferromag- 
netic Hund's rule coupling on the FCS for the orbital- 
degenerate Anderson impurity in the particle-hole sym- 
metric case. Using the RPT, we have derived the CGF 
for nonequilibrium-current distribution. The result is 
asymptotically exact at low energies and is described by 
the quasiparticles of the local Fermi liquid. Specifically, 
the explicit expression of the Fano factor for the shot 
noise is obtained in the fully screened higher-spin Kondo 
limit, which depends only on the orbital degeneracy M . 
We have also investigated the crossover between the large 
and small Hund's coupling limits in the two-orbital case 
(M = 2), using the NRG approach. Furthermore, the 
CGF indicates that the Hund's coupling gives rise to sin- 
glet and triplet pairs of quasiparticles carrying charge 2e 
in the backscattering current and these correlated charges 
characterize the current fluctuation through the dot. 
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